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Abstract 

In the last few years, lattice QCD has made a dramatic progress in understanding 
the physics of hadrons containing heavy quarks, from the first principle. This review 
summarises the maior achievements. 
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On . i. INTRODUCTION 

. Heavy quarks play a dominant role in the understanding of the weak interactions as 

well as of what may lie beyond the standard model. Since quarks are always bound into 
hadrons, and we do not understand the strong interactions rigorously, the results of these 
phenomena are often expressed in terms of non-perturbative parameters (also called matrix 
elements) that reflect our knowledge/ignorance of the strong interaction effects. Lattice 
QCD offers the best route to a non-perturbative determination of these parameters. 
Qh! The field of lattice gauge theories has nowadays achieved a maturity level where 

statistical errors are beaten down enough to expose systematic effects. This control over 
systematic effects has helped reliable extraction of matrix elements that allow us to test the 
Standard Model. In fact the results of these numerical experiments are now acknowledged 
by the Particle Data Group Q. What I present below are excerpts from results of many 
lattice QCD calculations; more details can be found in the proceedings of the annual lattice 
field theory conferences || ||. 

Lattice QCD simulations of hadrons containing heavy quarks provide a testing ground 
for many new ideas. The heavy hadrons are smaller than the light hadrons, and hence less 
susceptible to effects of finite lattice volume and dynamical quarks. (Typical lattice size in 
simulations nowadays is 1.5 — 2 fm.) Their nonrelativistic nature allows an easy comparison 
of the results with those from phenomenological quark potential models. Moreover, the 
results do not suffer from the uncertainties of chiral extrapolations. 

Three different formalisms for heavy quarks have been used in lattice QCD to get at 
the properties of heavy hadrons: (1) static, (2) non-relativistic (NRQCD), (3) conventional. 



L x = Qv„D»Q (1) 

L 2 = Q\D° tjD 2 + -^-a ■ B + 0(^—)]Q (2) 
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L 3 = q(p+M)q (3) 



Here Q refers to 2— component heavy quark fields, defined so as to remove the static mass 
M from the Lagrangian (otherwise M would appear as an additive constant). The field q 
contains both particle and anti-particle components, and M cannot be removed in that case. 
The coefficient c\ has to be non-perturbatively determined by a matching between NRQCD 
and conventional formalisms. 

In the static case, the 4— velocity of the heavy quark is a constant of motion, and the 
calculations are often carried out in the rest frame. The formalism is meaningful when the 
hadron has only one heavy quark, and the results provide the limiting behaviour as M — > oo. 
The NRQCD formalism is organised in powers of M~ l for hadrons having only one heavy 
quark, and in powers of v ~ p/M ~ a s (p) for hadrons having more heavy quarks Q. 

The lattice momentum cutoff in typical simulations nowadays is n/a ~ 5 — 10 GeV. 
The NRQCD approach is useful for Ma > 1. With the inclusion of 0(M~ 1 ) terms, L2 
works well for b— quarks and can be extrapolated to the case of c— quark. The conventional 
approach requires Ma < 1. It works reasonably for c— quarks, and can be extrapolated to 
heavier masses. 

A good consistency check on the methods is that the NRQCD and the conventional re- 
sults agree in the region between mj and m c . This comparison requires various M— dependent 
renormalisations to justify the extrapolations in the quark mass, but at the end we have con- 
tinuous functional behaviour for the quantities measured. 

There are many technical details involved in converting lattice results to numbers 
useful for phenomenology. The techniques are sophisticated, and the caveats should be kept 
in mind, instead of blindly using the numbers in phenomenological analyses. I only mention 
these without elaboration: 

• Both the action and the operators are improved, using higher order terms in a and 
M _1 to cut down discretisation errors. The creation and annihilation operators for the 
asymptotic states are smeared to increase overlap with the desired wavef unctions. 

• The lattice perturbation theory is not well-behaved in terms of the bare coupling. 
Rather one needs to express the results in terms of the renormalised coupling that is 
extracted from the simulation itself. A poor man's choice is "tadpole-improved" per- 
turbation theory. With a non-perturbative cutoff, the lattice results are unambiguous, 
but renormalon ambiguities appear when converting them to a continuum scheme such 
as ~MS. 

• The symmetries of the lattice theory are only a subset of the continuum ones. This 
makes renormalisation and mixing of the lattice operators more complicated than in 
the continuum. In principle, all the renormalisation constants should be determined 
through non-perturbative definitions. In practice, it turns out that finite and loga- 
rithmic renormalisations (e.g. dim< 4 operators such as quark bilinears) can be well 
estimated using the renormalised perturbation theory. Power law divergences (e.g. 
dim> 4 operators mixing with lower dimension operators), however, have to be elimi- 
nated using non-perturbative subtractions based on Ward identities. 



• A majority of lattice results, and certainly the ones where all the errors are well under 
control, are obtained in the quenched approximation ignoring the effects of dynamical 
quarks. This approximation violates the unitarity of the theory and unphysical sin- 
gularities appear in the chiral limit. Simulation results indicate that these violations, 
however, are not serious for typical simulations carried out with the light quark about 
as heavy as the strange quark. For certain quantities, results including two or four 
flavours of dynamical quarks have become available. These are combined with the 
quenched (i.e. rif = 0) values to obtain physically relevant numbers (i.e. n/ = 3). 

• Specific scaling forms are used to extrapolate the results to the infinite volume, con- 
tinuum and chiral limits. Algorithms for simulation of the path integral and for cal- 
culation of the quark propagator in a background gluon field keep on getting revised. 
Parametrisation of the highly correlated lattice data needs careful fitting procedures. 

II. PROPERTIES OF HEAVY-HEAVY MESONS 

Spectroscopy of heavy quarkonia (Upsilon, Charmonium and B c ) is the simplest test 
of lattice simulations. There are two bare parameters: the coupling and the heavy quark 
mass. Effects of light quarks appear only through vacuum polarisation loops. NRQCD 
simulations directly extract the level splittings, and the lattice cutoff is determined (in GeV) 
by assigning one of the splittings its physical value. The splittings are largely insensitive to 
the precise value of M, which keeps systematic errors under control. A remarkable feature of 
the results is that many of them have been obtained on rather coarse lattices using improved 
actions. 

Fig.l shows how well the lowest spin- averaged states of the Upsilon spectrum can 
be reproduced ||. A closer look reveals that the effects of light dynamical quarks shift the 
numbers closer to their experimental values. The Upsilon hyperfine splittings, shown in Fig.2, 
are more sensitive to the details of the lattice action: quenching tends to underestimate them, 
and accurate values are needed for M and C\ (which undergoes substantial renormalisation). 

Results for Charmonium and B c (not yet discovered) spectra are not as impressive. 
The main difficulty is that the quarks are more relativistic and higher order terms in the 
NRQCD expansion are needed to keep the systematic errors under control. The results 
are still phenomenologically useful and can be used to tune quark models. An example of 
Coulomb gauge wavefunctions for Charmonium states is shown in Fig. 3 ||. 

Calculations of heavy-light hadron masses suffer from the additional uncertainty of 
extrapolating the light quark mass to the chiral limit. Results are again best expressed as 
level splittings. Reasonable values are found for B s — B, At, — B and £& — A& splittings. 
Hyperfine splittings (e.g. B* — B and — £ 6 ), however, are typically underestimated and 
require more care. 

First results for annihilation decays of quarkonia have recently become available [[?[]. 
Annihilation decay rates, T(QQ — > light hadrons, , yj,l + l~), can be expressed as a power 
series in M _1 , with each term factorised into a short distance perturbative coefficient and a 
non-perturbative hadronic matrix element of 4— fermion operators ||. In the vacuum satu- 
ration approximation (VSA), the "colour singlet" combinations can be written in terms of 



quarkonium wave functions at the origin and their derivatives. The "colour octet" contri- 
butions measure the QQg component in the P— wave states. The numerical results are in 
reasonable accord with the VSA and the experimental data. 

III. THE QCD RUNNING COUPLING 

To verify that QCD is indeed the theory of the strong interactions, it must be demon- 
strated that perturbative and non-perturbative determinations of the coupling a s agree. The 
former is typically extracted from experimental results of high energy processes. The latter 
can be evaluated on the lattice from low energy hadronic properties. To this end, several 
non-perturbative definitions of ot s have been employed (note that otjjg does not exist outside 
perturbation theory): 



a q - q (q =-) = C^r 2 F(r) ^ oo (4) 
r 

o<-sf{<1=j) oc [-Q-^Z(L,r])] v= o ^ (5) 

1 (P x (0)P x (L)) g ^o + . ... 

aMq = l ] x (p z (o)p s {l)) constant (6) 

f 3.4 3 ln(^i x i) (7 . 

ap{q= - ] = "4^(1 -(1.19 + 0.07^)^ °° (7) 

The infrared behaviour of these definitions differ, as indicated above, but they should all come 
together towards the ultraviolet end. Strictly speaking, with a non-perturbative regulator, 
scaling is exact only on the critical surface. Away from g — 0, scaling holds only upto a 
certain tolerance level, and different definitions of a s help quantify this tolerance level. A 
comparison of g 2 SF and g\ P for the pure gauge SU{2) theory, shown in Fig. 4 ||, gives us 
faith to use perturbative scaling in the actual region of simulations. 

Having demonstrated scaling, the lattice results can be converted into a value for 
ajfg(Mz)- The common procedure is to use the spectrum results for quarkonia to determine 
ap at a specific value of the cutoff, and then use perturbation theory to convert it into 
a JTs- Given that a s = 0(0.1), a precision of 0(1%) requires 2— loop perturbative matching. 
Determinations of a s from different spectral quantities need not agree with each-other, when 
the number of dynamical quarks is unphysical. It is reassuring to see, as illustrated in Fig. 5 



T0[| , that they indeed converge when extrapolated to n/ = 3. Results obtained by various 
groups and their extrapolations to n/ = 3 are displayed in Fig.6 ||; they typically lie at the 
lower edge of the band of values from non-lattice determinations. 

There is yet another step in converting lattice numbers to a physical a-jjg(Mz)- The 
dynamical quarks included in lattice simulations are much heavier than the physical u and d 
quarks, and do not correctly reproduce the effects of pion loops. Chiral symmetry breaking 
in QCD implies that extrapolation of spectral quantities to the chiral limit be linear in the 
quark mass flli"|, in contrast to quadratic dependence seen in weak coupling perturbation 



theory. The final result is 



a m (M z ) = 0.117(1)(1)(2) 



(8) 



where the first error combines statistical and discretisation uncertainties, the second is due 
to vn^yni and the third comes from perturbative conversions. 



IV. PROPERTIES OF HEAVY-LIGHT PSEUDOSCALARS 

The pseudoscalar decay constant is a phenomenologically important parameter. Though 
its calculation on the lattice is straightforward, the number has remained hard to pin down 
due to its substantial mass dependence — Fig.7 [12| shows the extent of 0{M~ V ) corrections 



to the static limit. Exploiting the freedom that the various 0(M~ l ) contributions can be 
selectively extracted from lattice simulations, the high sensitivity of /q 9 to quark mass has 
been traced down to the large heavy quark kinetic energy |T3 |. Fig. 8 |T2j shows the ex- 



trapolation of f B values to the continuum, and also indicates that dynamical quarks would 
increase f B somewhat. Overall it is safe to quote O] 

f B = 180(20) MeV (with f n = 132 MeV) . (9) 

The ratios of decay constants are more stable as several correlated systematic errors cancel, 
and the results are [fT2 , |i~4 |: 



f B /f D « 0.85 , f B Jf B « 1.10 , fnjfj, « 1.10 , (10) 

all with errors of about 5%. 

Lattice calculations have also estimated several bound state parameters often used in 
B— physics phenomenology. In absence of any asymptotic quark states in QCD, these are 
defined through behaviour of lattice quark propagators in the Landau gauge. To quote [|T^] : 



the B— meson binding energy A = 180(30) MeV, the MS running quark mass m^mj,) = 
4.15(20) GeV, the heavy quark kinetic energy in B— meson Ai = 0.1(1) GeV 2 . 

The B— parameter characterising B — B mixing is another useful quantity. In the 
effective theory below mj, it has no anomalous dimension at 1— loop. It can be easily seen, 
at least at strong coupling, that the deviation of B B from one is 0(M~ 1 , N~ 2 ). Indeed 
different lattice calculations confirm that B B (mb) = 0.90(5) [[UJ. Dynamical quarks do not 
seem to have substantial effect on B B . Increasing the light quark mass increases B B slightly, 
so that Bb s > Bbd- Conversion of the lattice B B to the RG-invariant value used in continuum 
analyses is rather uncertain: B B = 1.3(2). 

V. FORM FACTORS OF HEAVY-LIGHT HADRONS 

The spin-flavour symmetries of heavy quarks reduce form factors of heavy-light hadrons 
to Isgur-Wise functions |T7|]. Matrix elements of the charge operators (i.e. pi = pf) deter- 



mine the form factor normalisations, and provide scaling rules for their M— dependence. In 
the heavy quark effective theory, the charge operators act in the rest frame, and hence con- 
strain the form factors at q 2 = q^^- On the other hand, the charge operators correspond to 
q 2 = in the lightcone frame, and thus the lightcone sum rules constrain the form factors 
at maximum recoil. The two complementary constraints can be merged together in pole 
dominance ansatze for the form factors, where the crucial property is the order of the pole. 
Lattice calculations can be used to verify these ansatze. In practice, it is easier to measure 



form factors near g^ ax in experiments, and near q 2 = on the lattice. The pole dominance 
ansatze are often conveniently adopted to extend them to other q 2 . 

Lattice calculations can observe how the symmetry relations amongst the form fac- 
tors get modified as the quark masses are reduced. For the semileptonic pseudoscalar-to- 
pseudoscalar decay, the pole relation between /+ and fo is well satisfied, even when the 
the heavy quark mass is 0(m c ), as exemplified in Fig.9 |18| . The (9(M _1 ) corrections to 



the relations amongst the form factors appearing in the semileptonic pseudoscalar-to-vector 
decay are larger, in particular the relation between V and Ai works poorly even at m^. The 
corrections seem to affect mostly the normalisations, however, and not the q 2 — dependence. 
The form factors occuring in the radiative decay B — > K*j, and heavy-to-light semileptonic 
decays B — > plv and B — > nlu, can be related using heavy quark and light flavour SU(3) 
symmetries. These relations seem to work well, but the q 2 — dependence of the form factors 



is not yet fully sorted out in lattice calculations [16, 19 



The Isgur-Wise functions have been extracted from lattice form factors evaluated at 



sufficiently heavy quark mass. Examples for meson and baryon cases are shown in Fig. 10 []I 
and Fig.ll |2D| respectively. The meson Isgur-Wise function, £(t>j • v /), shows an increase as 
the light quark mass is reduced, and the slope at zero recoil is £'(1) ~ —1 in the chiral limit. 
First results for the baryon Isgur-Wise function |2(J indicate that it falls off more steeply 



than the meson one, and its slope decreases by about a factor of two between m s and the 
chiral limit. 

In strong coupling lattice QCD, £(t>; • Vf) has a double pole form, with the singularity 
at Vi ■ Vf = — 1 This result can be rephrased as PCVC: the divergence of the partially 



conserved vector current is saturated by the scalar meson pole near Mi + Mf (the scalar pole 
may arise from the combination of a set of states spread over an energy range much smaller 
than M). This ansatz fits both the lattice and experimental (Fig. 12 P3 shows the latest 



results) data as well as any other fit. It is worth exploring if there is any limit of QCD, quite 
likely it would involve N c — > oo, where this ansatz becomes exact. 

VI. CONSTRAINING THE STANDARD MODEL 

Lattice results can help in precision tests of the standard model, by pinning down the 
elements of the Cabbibo-Kobayashi-Masakwa quark mixing matrix. Of particular interest is 
the unitarity triangle involving the third generation, which determines CP— violation. The 
present constraints can best be expressed in terms of the Wolfenstein parameters p and rj, 
and they are displayed in Fig. 13. Particularly useful here are the lattice inputs for Bb 
and Bk [p3| . The bounds suggest that, within the standard model, CP— violation is close 
to maximal (i.e. the phase 5 ~ tt/2). The standard model will be constrained further by 
results leading to |e'/e|. 

Lattice B B values also determine, in terms of |Vt s |, the extent of B s — B s oscilla- 
tions. These oscillations may be observable at LEP, SLC and HERA-B |19| . 
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Figure 1: T spin- averaged spectrum [5]. 
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Figure 2: T spin-splittings [5]. The symbols have the same meaning as in Fig.l. 
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Figure 3: The radial wavefunctions for the IS and IP Charmonium states computed using 
improved actions and two different lattice spacings [6]. Wavefunctions from a continuum 
quark model are also shown. Statistical errors are negligible for the IS wavefunction. 
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Figure 4: A comparison of the pure gauge SU(2) running coupling for two different 
non-perturbative definitions [9]: squares are gg F and circles £1X6 Q'j'p- 

The energy scales 

have been matched with the appropriate ratio of the A— parameters. The dotted line shows 
the perturbative 2— loop evolution. 
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Figure 5: a"p f \s.2GeV] versus rif [10]. The scale is set by the T 1S-1P (full circles), the T 
1S-2S (open circles) or the Charmonium 1S-1P (boxes) splitting. 
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Figure 6: a P ' versus ln(^) [5]. The symbols are: rif — charm (□), rif — bottom (O), 
rif — 2 charm (+), rif — 2 bottom (x), and rif — 3 extrapolated (fancy boxes). The region 
between the two full curves corresponds to aj^(Mz) G [0.115,0.125], typically obtained in 
non-lattice determinations. The dotted lines are 3— loop scaling curves for rif — 0. 
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Figure 7: fQ q (MQ q )2 vs. 1/Mq q for rif = 0, [3 = 6.3 Wilson fermions [12]. The solid line is 
a quadratic fit to the octagons ("heavier heavies" + static); the dashed line is a quadratic 
fit to the crosses ("lighter heavies" + static). 
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Figure 8: vs. a for nj = (octagons) and rif = 2 (crosses) data [12]. The solid line is 
a linear fit to all rif — points; the dashed line is a constant fit to the three rif — points 
with a < 0.5 GeV" 1 . The extrapolated values at a = are indicated by bursts. The scale is 
set by f n = 132 MeV throughout. 
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Figure 9: Comparison of semileptonic decay form factors for the heavy-light pseudoscalar 
meson [18]. The heavy quark mass roughly changes from m c to m s , while the light quark 
mass is about 0Am s . 
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Figure 10: £(w)/£(l) at various quark masses [18]. The initial heavy quark mass is about 
m c , and various symbols label the masses it decays into (from m c to m s ). The clusters of 
three points show an increase in the form factor with a decrease in the spectator quark mass 
(from m s to 0Am s ). 




Figure 11: The baryon Isgur-Wise function from the elastic form factor [20]. Different 
symbols denote different heavy quark massess in the vicinity of m c , while the light quark 
masses straddle m s . The points with the same u depict an increase in the form factor with 
a decrease in the spectator quark masses. 
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Figure 12: The B — > D*W form factor from CLEO II [22]. y — V{ • Vf and F(y) reduces to 
the Isgur-Wise function in the limit m&,m c — > oo. Combination of ARGUS, CLEO II and 
ALEPH results, with the latest world average lifetimes [1], gives the linear fit parameters: 
V cb F(l) = 0.0359(22) and slope p 2 = 0.88(12). 




Figure 13: Constraints on the Wolfenstein parameters of the quark mixing matrix, from the 
combination of current experimental [1] and lattice results. The solid semi-circles are from 
|Kt&/Vd,|, the dashed hyperbolas are from |e| and Bk, and the dotted circles are from Arris 
and f%B B . 



